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I.  INTRODUCTION 


Surface  irregularities  such  as  cavities,  bumps  and  flares,  etc.,  may  all 
be  present  in  actual  missile  configurations.  The  resulting  flow  field  for 
these  shapes  Is  characterized  by  locally  separated  flow  regions  due  to  these 
surface  irregularities.  Additionally,  large  separated  flow  regions  exist 
behind  the  base  of  the  missile.  These  separated  flow  regions  could  have  a 
large  affect  on  the  missile  aerodynamics.  Fins  and  control  surfaces  are  often 
added  to  Improve  the  stability.  Proper  location  of  these  devices  are  crucial 
in  achieving  full  effectiveness.  The  control  effectiveness  of  the  fins  or 
other  control  surfaces  can  be  seriously  degraded  if  they  are  located  in 
regions  of  large  separated  flow.  A  knowledge  of  the  expected  flow  field  and 
the  boundary  layer  properties  can  thus  be  very  important  in  the  initial  design 
phase.  The  Aerodynamics  Technology  Branch  of  the  US  Army  Missile  Laboratory 
has  requested  numerical  flow  field  data  for  a  new  missile  design.  Numerical 
computations  have  been  made  for  the  requested  missile  geometry  at  M  =  0.6. 
The  results  of  these  computations  are  presented  in  this  memorandum  report. 


II.  COMPUTATIONAL  TECHNIQUE 

The  Azimuthal  Invariant  (or  Generalized  Axi symmetri c)  thin-layer  Navier* 
Stokes  equations  for  general  spatial  coordinates  5,  n,  c  can  be  written  as1 


-I. 
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where  £  *  E(x,y,z,t)  is  the  longitudinal  coordinate 
n  *  n(y,z,t)  is  the  circumferential  coordinate 
C  3  c(x,y,z,t)  is  the  near  normal  coordinate 
t  3  t  Is  the  time 
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the  Azimuthal- Invariant  Navien-Stok.ee  Equations,"  U.S.  Army  Ballistic 
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'■torch  1980.  (AD  A08S7161  (Also  3ee  AIAA  Paper  No.  79-0010 ,  ■ January 
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The  velocities 

U  *  ?t  +  5xu  +  V  +  5z" 

V  =  nt  +  nxu  +  nyv  +  n2«  (2) 

W  *  Ct  >  txu  +  Cyv  +  c2« 

represent  the  contravariant  velocity  components. 

The  Cartesian  velocity  components  (u,  v,  w)  are  nondimensionalized  with 
respect  to  am  (free  stream  speed  of  sound).  The  density  (p)  is  referenced  to 
t>m  and  total” energy  (e)  to  p^a2.  The  local  pressure  is  determined  using  the 
equation  of  state. 


P  »  (y  -  l)[e  -  Q.5p(u2  +■  v2  +  w2)] 


where  y  is  the  ratio  of  specific  heats. 

In  Equation  (1),  axisymmetric  flow  assumptions  have  been  made  which  result 

in  the  source  term,  H.  The  details  of  how  this  is  obtained  can  be  found  in 
Reference  1  and  are  not  discussed  here.  Equation  (I)  contains  only  two 
soatial  derivatives.  However,  it  retains  all  three  momentum  equations  and 
allows  a  degree  of  generality  over  the  standard  axisymmetric  equations.  In 


particular,  the  circumferential  velocity  is  not  assumed  to  be  zero,  thus 
allowing  computations  for  spinning  projectiles  to  be  accomplished. 


The  numerical  algorithm  used  is  the  Beam-Warming  fully  implicit,  approxi¬ 
mately  factored  finite  difference  scheme.  The  algorithm  can  be  first  or 
second  order  accurate  in  time  and  second  or  fourth  order  accurate  in  space. 
Since  the  interest  is  only  in  the  steady-state  solution,  Equation  (1)  is 
solved  in  a  time  asympotic  fashion  and  first  order  accurate  time  differencing 
is  used.  The  spatial  accuracy  is  fourth  order.  Details  of  the  algorithm  are 
included  in  References  2-4. 

For  the  computation  of  turbulent  flows,  a  turbulence  model  must  be  sup¬ 
plied.  In  the  present  calculations  a  two  layer  algebraic  eddy  viscosity  model 
by  8aldwin  and  Lomax5  is  used.  In  their  two  layer  model  the  inner  region 
follows  the  Prandtl-Van  Oriest  formulation.  Their  outer  formulation  can  be 
used  in  wakes  as  well  as  in  attached  and  separated  boundary  layers.  In  both 
the  inner  and  outer  formulations,  the  distribution  of  vorticity  is  used  to 
determine  length  scales,  thereby  avoiding  the  necessity  of  finding  the  outer 
edge  of  the  boundary  layer  (or  wake).  The  magnitude  of  the  local  vorticity 
for  the  axi symmetric  formulation  is  given  by 
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In  determining  the  outer  length  scale  a  function5 


F(y)  =  yM  [I  -  exp(-y+/A+)] 


(5) 


was  used. 

The  thin-layer  Navier-Stokes  computational  technique  described  above  was 
used  in  conjunction  with  an  unique  flow  field  segmentation  nrocedure  and  a 
computational  capability6  7  has  been  recently  developed  for  predicting  the 


2.  -7.  £.  Stager,  "Implicit  Finite  Difference  Simulation  of  Flow  About 
Arbitrary  Geometries  with  Application  to  Airfoi Is . "  .4I7L4  J ouma l ,  V ol.  IS, 
■Vo.  7,  July  1973,  pp.  679-696, 

3.  T.  H.  Pulliam  and  J.  L.  Steger,  "On  Implicit  Finite-Difference  Simulations 
of  Three-Dimensional  Flow,"  AIAA  Journal,  Vol,  Id,  So.  2,  February  1980  , 
pp.  159-167. 

4.  R.  3eam  and  R.  F.  Warning,  "An  Implicit  Factored  Scheme  for  the  Compres¬ 
sible  Savier-Stokes  Equations,"  AIAA  Journal,  Vol.  16,  So.  4,  April  1973, 
pp.  393-402. 

5.  3.  S.  Saldwin  and  H.  lomas,  'Thin-Layer  Approximation  ana  Algebraic  'lodei 
for  Separated  Turbulent  Flows,"  AIAA  Paper  So.  73-237 ,  1979. 
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flow  field  over  a  projectile  including  the  base  region  flow.  The  flow  field 
segmentation  procedure  allows  the  complete  numerical  simulation  of  a  projec¬ 
tile.  The  details  of  this  procedure  can  be  found  in  References  6  and  7.  This 
code  is  used  here  to  calculate  the  full  flow  field  over  a  missile  at  M  *  0.6 
and  o*0. 


III.  RESULTS 

Numerical  computations  have  been  made  at  two  Reynolds  Numbers,  Re  *  4.2 
x  106/ft  and  3.0  x  106/ft  for  turbulent  flow  and  Re  =  4.2  x  10 6/ft  for 
laminar  flow  (total  length,  L  *  23.63").  Solutions  were  marched  in  time  until 
the  steady  state  results  were  achieved.  The  lengths  shown  in  results  are  in 
calibers  (1  caliber  3  2.66",  the  reference  diameter). 

Figure  1  shows  an  expanded  view  of  the  computational  grid  near  the  missile 
which  is  approximately  9  calibers  long.  The  grid  consists  of  201  points  in 
longitudinal  direction  and  50  points  in  the  normal  direction.  The  dark 
regions  result  from  grid  clustering.  These  are  the  regions  where  the  flow 
variables  are  expected  to  change  considerably.  The  expanded  grid  in  the  base 
region  is  shown  in  Figure  2  and  has  been  adapted  to  the  wake  shear  layer  as 
the  solution  developed.  Computed  results  are  now  presented  for  the  various 
cases. 

A. 

Figure  3  shows  the  velocity  vectors  in  the  regions  where  surface  irregu¬ 
larities  are  present  and  flow  separation  may  occur.  As  seen  in  Figure  3(a) 
there  is  a  very  thin  region  of  reversed  flow  near  the  surface.  Figure  3(b) 
and  (c)  are  expanded  views  of  Figure  3(a).  Figure  3(d)  shows  the  velocity 
vectors  near  the  flare  and  again  only  a  very  thin  separated  flow  region  is 
predicted  near  the  missile  surface. 

Qualitative  features  of  the  base  region  flow  field  are  shown  in  Figures  4- 
7.  Figure  4  is  a  plot  of  the  velocity  vectors  and  clearly  shows  the  recircu- 
lator.y  flow  in  the  near  wake.  Figure  5  shows  the  particle  paths  in  the  base 
region  and  more  clearly  shows  the  separation  bubble.  Mach  number  and  normal¬ 
ized  static  temperature  contour  plots  are  shown  in  Figure  6  and  7,  respec¬ 
tively,  and  show  the  strong  free  shear  layer.  It  is  clear  from  Figure  6  that 
the  flow  is  subsonic  everywhere  in  the  base  region. 


S.  J .  Saku,  C.  J.  Hietubiaz  and  J.  L.  Steger,  ’'Numerical  Computation  of  3ase 
Flow  for  a  Projectile  at  Transonic  Speeds ,”  U .3.  Army  Sallistic  Research 
Laboratory ,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-02495 ,  June  1993. 
(AD  A130293)  (Also  see  AIAA  Paper  Ho.  92-1338,  August  1982.) 

7.  J.  Saku,  C.  J.  Hietubiaz  and  J.  L.  Steger,  "Havier-Stokes  Computations  of 
Projectile  Base  Flou  with  and  without  Base  Injection,”  U.3.  Army  Ballistic 
Research  Laboratory ,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-02532 , 
November  1983.  (AD  A13S738)  ( AI30  see  AIAA  Paper  Ho.  33-0224,  January 
1983.) 


Case  1 ,  Re  »  4.2  x  IQ6/ ft  (Turbulent  Flow) 


Quantitative  results  are  presented  in  the  form  of  pressure  coefficient  and 
skin  friction  coefficient.  Figure  8  shows  the  surface  pressure  distribution 
as  a  function  of  the  longitudinal  position.  The  flow  is  subsonic  everywhere 
(Cp  >  Cp*)  and  the  expansions  and  recompressions  are  observed  at  the  cavity, 

protuberance  and  flare  as  expected.  The  pressure  distribution  shown  down¬ 
stream  of  the  base  corner  is  along  a  line  followng  the  shear  layer.  Figure  9 
shows  the  skin  friction  coefficient  distribution  and  clearly  indicates  the 
reversed  flow  regions  where  Cf  <  0. 

B. 

The  results  for  this  case  are  essentially  the  same  as  those  of  the 
previous  case.  Small  differences  in  the  surface  pressure  distribution  and 
skin  friction  coefficient  distribution  are  observed  and  these  distributions 
are  shown  in  Figures  10  and  il.  The  differences  are  very  small  and  usually 
occur  in  the  separated  flow  regions. 

C. 

Qualitative  features  of  the  base  region  flow  field  for  the  laminar  flow 
case  are  shown  in  Figures  12-15.  Although  small  differences  can  be  observed, 
the  extent  and  size  of  the  wake  separation  bubble  remains  essentially  the  same 
as  for  the  turbulent  case.  Significant  differences  in  the  flow  field,  how¬ 
ever,  occur  in  the  regions  of  the  cavity,  protuberance  and  the  flare  and  are 
shown  in  Figures  16-19.  Figures  16  and  17  show  the  velocity  vectors  in  these 
regions  and  the  large  separated  flow  regions  are  evident.  Figures  18  and  19 
are  stream  function  contour  plots  for  the  same  regions  and  more  clearly  show 
the  reversed  flow.  These  stream  function  plots  clearly  show  the  separation 
and  reattachment  points. 

Quantitatively,  the  surface  pressure  and  the  skin  friction  coefficient 
distributions  are  shown  in  Figures  20  and  21,  respectively.  As  seen  in  Figure 
20,  the  expansions  and  the  compressions  are  weak  compared  to  the  turbulent 
case.  The  magnitude  of  the  skin  coefficient  shown  in  Figure  21  is  consider¬ 
ably  lower  than  for  the  turbulent  case  (Figure  9).  Additionally,  larger 
regions  of  separated  flow  are  clearly  seen  (Cf  <  0). 

As  part  of  the  missile  design,  fins  are  located  at  X/D  *  8.306.  The  flow 
field  at  this  location  is  of  considerable  importance  to  the  Aerodynamic 
Technology  3ranch  of  the  US  Army  Missile  Command.  The  flow  field  requested  is 
in  terms  of  plots  and  tabular  data  of  density,  velocity,  static  pressure  and 
dynamic  pressure  profiles.  Figures  22-26  show  the  profiles  of  density,  velo¬ 
city  components  in  streamwise  and  normal  directions,  static  and  dynamic  pres¬ 
sures  respecti vely.  The  tabular  data  of  these  profiles  are  included  in  Table 
1  and  includes  the  following  variables: 


Z/D  normal  distance  measured  from  the  surface 

RH0/RINF  non-dimensional  density 

U/UINF  non-dimensional  x-component  of  velocity 

W/U INF  non-dimensional  y-component  of  velocity 

P/PINF  non-dimensional  static  pressure 

Q/QINF  non-dimensional  dynamic  pressure 


Case  2,  Re  3  3.0  x  106/ft  (Turbulent  Fit 


where  0  Is  the  reference  diameter  and  RINF,  UN IF,  PINF  and  QINF  are  the  free 
stream  values  for  density,  velocity,  static  pressure  and  dynamic  pressure 
respectively. 


IV.  CONCLUDING  REMARKS 

A  thin-layer  Navi er- Stokes  code  has  been_  used  to  compute  the  entire  flow 
field  about  a  missile  configuration  of  interest  to  the  Aerodynamics  Technology 
8ranch  of  the  US  Army  Missile  Laboratory.  Numerical  computations  have  been 
made  at  Mw  *  0.6  and  a  *  0  and  for  both  laminar  and  turbulent  flow  condi¬ 
tions.  Turbulent  flow  computations  indicate  very  thin  regions  of  separated 
flow  near  the  missile  surface  in  the  neighborhood  of  the  Irregularities  in 
surface  geometry.  The  laminar  case  predicts  much  larger  separated  flow  re¬ 
gions  In  the  same  areas.  The  flow  field  in  terms  of  plots  and  tabular  data  of 
density,  velocity,  static  pressure  and  dynamic  pressure  at  the  location  of 
fins  are  also  included. 
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Figure  2.  Expanded  View  of  the  Grid  in  the  Base  Region 
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Figure  17.  Velocity  Vectors,  M„  *  0.6,  a  ■  0,  Re  -  4.2  x  106/ft  (Laminar) 
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Table  1.  Flow  Field 


Z/D 

RHO/RIHF 

O/tJIHF 

0.00000 

1.03619 

0.00000 

.00004 

1.03321 

.01940 

.00008 

1 .02994 

.04105 

.00014 

1.02614 

.06646 

.00022 

1.02179 

.09581 

.00031 

1.01705 

.12779 

.00044 

1.01223 

.16018 

.00060 

1.00763 

.19064 

.00080 

1.00347 

.21754 

.0010T 

.99982 

.24027 

.00141 

.99671 

.25912 

.00186 

.99411 

.27488 

.00243 

.99203 

.28850 

.00316 

.99043 

.30083 

.00412 

.98925 

.31256 

.00534 

.98842 

.32427 

.00692 

.98787 

•33641 

.00896 

.98757 

.34938 

.01159 

.98750 

.36353 

.01498 

.98763 

.37912 

.01935 

-98791 

.39634 

.02498 

.98831 

.41548 

.03225 

.98883 

.43716 

.04160 

.98952 

.46264 

.05365 

.99052 

.49394 

.06918 

.99203 

.53343 

.08919 

.99419 

.58291 

.11505 

.99705 

.64250 

.14870 

1 .00043 

.70988 

.19295 

1 .00391 

.78045 

.25183 

1 .00689 

.84828 

.33023 

1.00881 

.90762 

.43323 

1.00928 

.95441 

.56605 

1.00828 

.98713 

.73528 

1.00618 

1 .00682 

.95004 

1 .00359 

1.01638 

1 .22291 

1.00118 

1.01942 

1 .57059 

-99939 

1.01919 

2.01496 

.99842 

1.01790 

2.58437 

.99817 

1.01668 

3-31445 

.99840 

1.01579 

4.24764 

.99881 

1.01509 

5.43206 

.99917 

1.01431 

6.92126 

.99938 

1.01325 

8.77616 

.99943 

1.01181 

11.06881 

.99939 

1 .00998 

13.88611 

.99937 

1.00778 

17.33081 

.99946 

1.00531 

21.52252 

.99968 

1.00267 

26.60690 

1.00000 

1 .00000 

Oata  at  X/D  *  8.306 


W/UINF 

P/PMF 

Q/QINF 

0.00000 

1.03619 

0.00000 

•00399 

1.03619 

.00041 

.00845 

1.03619 

.00181 

.01368 

1.03618 

.00472 

.01971 

1.03617 

.00978 

.02628 

1.03615 

.01731 

.03291 

1.03612 

.02707 

.03911 

1.03605 

.03816 

.04450 

1.03594 

.04948 

.04893 

1.03577 

.06011 

.05247 

1.03557 

.06966 

.05533 

1.03538 

.07816 

.05777 

1.03528 

.08588 

.05996 

1.03532 

.09319 

.06196 

1.03543 

.10045 

.06376 

1.03556 

.10795 

.06532 

1.03565 

.11601 

.06663 

1.03571 

.12493 

.06768 

1.035T3 

.13503 

.06841 

1.03572 

.14657 

.06872 

1.03566 

.15986 

.06848 

1.03553 

.17524 

.06759 

1 .03532 

.19349 

.06604 

1.03504 

.21611 

.06390 

1.03469 

.24571 

.06131 

1 .03425 

.28601 

.05838 

1 .03358 

.34120 

.05515 

1.03246 

.41462 

.05158 

1 .03055 

.50681 

.04759 

1.02760 

.61376 

.04315 

1.02356 

.72641 

.03826 

1.01865 

.83251 

.03294 

1.01331 

.92045 

.02722 

1.00806 

.98324 

.02126 

1 .00334 

1.02041 

•01532 

.99948 

1.03698 

.00978 

.99667 

1.04054 

.00504 

.99493 

1.03813 

.00139 

.99415 

1.03449 

-.00111 

.99411 

1.03175 

-.00260 

.99452 

1.03019 

-.00336 

.99512 

1.02918 

-.00368 

.99569 

1.02798 

-.00375 

.99615 

1.02605 

-.00366 

.99653 

1.02319 

-.00340 

.99689 

1.01945 

-.00291 

.99738 

1.01500 

-.00214 

.99807 

1.01011 

-.00113 

.99897 

1.00503 

0.00000 

1.00000 

1.00000 
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